HYDRODYNAMIC TURBULENCE AS A PROBLEM 
IN NONEQUILIBRIUM STATISTICAL MECHANICS. 



by David Ruellef. 



^ ! Abstract: We reformulate the problem of hydro dynamic turbulence as a heat flow 

I problem. We obtain thus a prediction 

J C. = f-^inr(f + i) 

Q . 3 In K 3 

for the exponents of the structure functions ((|A^f |p) = r'^p). The meaning of the ad- 
justable parameter k is that when an eddy of size r has decayed to eddies of size r/n 
their energies have a thermal distribution. The above formula, with (In^)"^ = .32 ± .01 
I is compatible with experimental data. This agreement lends supports to our physically 

^ ■ motivated picture of turbulence. 
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Hydrodynamic turbulence is known to be a chaotic phenomenon [1,2,3]. This means 
that the time evohition (/*) of a turbulent fluid system belongs to a much studied class 
of deterministic dynamics with sensitive dependence on initial conditions [4,5,6]. The 
statistical properties of turbulence are described by an ergodic invariant state p for (/*) 
and, since chaotic dynamical systems have (uncountably) many ergodic states, a choice 
has to be made. A physically reasonable choice is that of so-called SRB states (see [7,8] 
and the references given there). 

It is fair to say that the chaotic nature of turbulence has been largely ignored by the 
turbulence community, and that the choice of an ergodic state to describe the statisti- 
cal properties of turbulence has been made by ad hoc assumptions (closure assumptions, 
Gaussianity, multifractal structure). Indeed, the study of SRB or "physical" states for 
the Navier-Stokes dynamics appears at first impossibly difficult. Nevertheless we propose 
here an approach of this sort: we bypass the mathematical problems of SRB states by 
using our understanding of the physics of a specific dynamical system, namely that corre- 
sponding to heat conduction, as seen from the point of view of nonequilibrium statistical 
mechanics. Our approach will thus use basic physical ideas, and approximations, rather 
than ad hoc assumptions. We shall in this manner obtain a surprisingly coherent view of 
the fluctuations in turbulence (intermittency). 

We shall concern ourselves with incompressible fluids in 3 dimensions, described by 
the Navier-Stokes equation, but without paying too much attention to the specific form 
of the dissipative term. The fluid, with velocity field v, will be enclosed in a cube Co of 
side £o, which we may consider for simplicity to have periodic boundary conditions. We 
choose an integer n > 1 and divide Cq — Cqi into cubes Cni of side In — ^qX'~^, with 
i = 1, . . . , K^^, where n is a positive integer. Let (pni be the homothety mapping Cq to 
Cni- One can choose 2(k/^ — 1) real vector fields 11^ on with J Ua — 0, div 1/^ = 0, and 
such that if the velocity field v satisfies J v = 0, divv = 0, there is a unique representation 



with Cnia £ R" This means that v has a wavelet decomposition into components (roughly) 
localized in the cubes Cni- 

We think now of the standard physical situation where energy is put into the fiuid 
at large spatial wavelength (i.e., small n) and dissipated at small spatial wavelength (i.e., 
large n). Intermediate values of n correspond to the inertial range, where the time evo- 
lution should in some sense be Hamiltonian. Specifically, Arnold [9] has shown how an 
inviscid fiow could be interpreted as geodesic flow on the group of volume preserving dif- 
feomorphisms. The corresponding Hamiltonian is the kinetic energy of the velocity fleld. 

We may thus think of the time evolution for (the flnitely many) coefficients Cnia as 
Hamiltonian, with external forces acting at low and high n. This is related to the physical 
concept of eddies as dynamical structures localized in space. However, instead of a cascade 
of eddies of smaller and smaller size, we think of a system of coupled Hamiltonian systems 
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which we can label (n, i). If we assume that the different systems (n, i) are weakly coupled, 
wc can reinterpret the global dynamics as a heat flow from the small n, where energy is 
input, to large n where it is dissipated, i.e., rapidly carried away to structures of molecular 
size of the fluid. Note that the multifractal description of eddy cascades [10,11,12] ignores 
interactions between (n,i), (n, i') except when these eddies are created from a common 
(n — 1, j). This corresponds to saying that the lateral interaction between the systems 
{n,i), {n,i') is weak, but this assumption does not appear to be essential in our approach. 

There is no hope for an exact study of the dynamics of the coupled systems {n,i), 
but we can get a first approximation from the Kolmogorov scaling theory of homogeneous 
turbulence [13]. Since this theory gives unique answers, the problem of selecting an SRB 
state does not occur here. According to the Kolmogorov theory, the fluid velocity corre- 
sponding to Cni is Vni ~ (ein)^^^ where e is the mean dissipation per unit volume, and the 
kinetic energy corresponding to Cni is 

(we have put the fluid density equal to 1); the corresponding temperature is 

k2iK^-l) k 2(k3_i) ^ > 

where k is Boltzmann's constant. In view of the value of /c, we see that T„ is huge for 
small n so that the flow of heat from high temperature to low temperature agrees with 
the energy cascade from small I to large £ in the fluid. Notice that the heat resistance 
{Tn — Tn+i)/e is very large, which agrees with a weak coupling between the systems (n, i) 
for different values of n. 

We see the situation as follows: a heat flow interpretation of the energy cascade in 
homogeneous turbulence is possible, using scaling laws, but ignores fluctuations (inter- 
mittency). To understand fluctuations we have to study the fluctuations of the energy 
flow in the Hamiltonian system of the coupled {n,i). This is a problem of nonequilib- 
rium statistical mechanics, a problem known to be difficult [14,15,16]. In general one 
would need the systems (n, i) to be chaotic in some sense (this is physically reasonable for 
3-dimensional hydrodynamics) but the Anosov assumptions of [15,16] are unreasonably 
strong. In the present situation, a rigorous analysis appears quite out of reach at this 
time. An approximate study is however possible, and will give more specific results than 
earlier multifractal approaches [10,11,12], and there will be a physical justification rather 
than ad hoc assumptions. 

While we have at this time no detailed understanding of heat fiow from the point 
of view of rigorous statistical mechanics, we expect that Fourier's law should hold under 
normal conditions. This is however no great help since Kolmogorov's theory yields the 
precise temperature distribution (1). As to the "microscopic" fiuctuations, they are a diffi- 
cult problem in nonequilibrium [17,18], being different in nature from the well-understood 
equilibrium fluctuations. Here we shall use the assumption that the systems (n, i) have 
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weak mutual coupling to justify a Boltzmannian energy distribution for each Hamiltonian 
system (n, i). 

Because of the large temperature gradient, the flow of energy is overwhelmingly from 
the system (n, i) to the systems (n + 1, j). To study this energy flow we use the conservation 
of energy and scaling as in the multifractal approaches [11] to write 

\Vnif/in = \v{n+l)jf/^n+l Or \v(n+l)jf = l^nif (2) 

Note that \v\'^ is proportional to the kinetic energy with a weight 1/ time spent in 

a certain spatial frequency range. We interpret then (2) to mean that, given the energy 
Vni = \vni\^ in {'IT', i), the velocity v = V(^n+i)j is fluctuating with Boltzmannian distribution 

I f I ^ 

Therefore the energy V = V(^n+i)j has the normalized distribution 

where we write from now on Vn instead if Vni, etc. We view (3) as an approximate, but 
physically motivated relation, the validity of which will be discussed below. Note that if 
we replace Vnhy Vn = n^Vn we have that V = V^+i is distributed according to 

1 V ~ 

cxr)(— — )dV 
V V 

We now discuss the structure functions, i.e., the moments 

for positive integer p, and the exponents C,p such that 

~ 4- or Cp ln£, ~ HV^^^) = -n • | In « + HV^'^) 

We have here 

-Vl/Vb r r ~ p-Vr.-l/Vr,-2 r ^ p-Vu/Vr,-l 



= / dV^"-^ ... / dVn^ v^/' 

J Vq J J Vn-2 J Vn-1 

r dVn'^^^^p^ ■ = v:Li r d^e-H'" = v:LiT{i + 1) 

^0 Vn-l Jo "J 

SO that by induction we find 



In conclusion we have the (approximate) prediction 



Using either the heat propagation or the eddy cascade picture, we see that k should be 
chosen such that the initial V^-distribution concentrated on one value for (n, i) thermalizes 
to values of Vn+i for the systems (n + 1, j) distributed according to 

This requires k sufficiently large. However, if the value of k is too large, several different 
temperatures will be present among the systems (n + connected with {n,j), and the 
y^_i_i-distribution will not be Boltzmannian. The picture we have in mind is a situation 
in Cjii which depends on the spatial wavelength: at wavelength of the order of the size of 
supp(t/Q, o (f)~l) a single value of the kinetic energy is present; the distribution broadens 
as the wavelength diminishes and becomes a thermal distribution when it is divided by k; 
at smaller wavelengths there are several patches with different temperatures. Of course 
a rigorous justification of this picture is well beyond the power of current mathematical 
methods. We can only claim this: k should be such that when an eddy of size r has 
decayed to eddies of size r/n their energies have a thermal distribution, after which the 
process can start again. In the dissipative range the distribution of Vn should be cut off 
at large Vn- Numerically, one finds that (4) fits the experimental data [19] well ([20] less 
well) with 1/logK = .32 ± .01, i.e., n between 20 and 25. Note that (4) gives ^3 = 1- 

From a physical point of view, one can try the following interpretation: the change of 
behavior as one passes from large to small wavelength corresponds to what is observed at 
the onset of turbulence. We use here the physical fact that transport is much faster in the 
turbulent than in the laminar regime. Therefore, when sufficiently small scales are reached 
and we are in the turbulent regime, thermalization takes place. We may compute the 
length ratio k = £n/in+i in terms of the Reynolds number Rc for the onset of turbulence. 
If e is the energy dissipation per unit volume and u the kinematic viscosity, the Kolmogorov 
length is ?7 = (z/"^/e)^/^, so that ly = (e?7^)^/'^. The velocity corresponding to the length 
A is given in the turbulent regime by vx — (eA)^/^. Therefore if the onset of turbulence 
corresponds to A, we have 

A^ A^ 1/4^ .^3/4^^3/4 

The critical Reynolds number i?c is not defined with precision, but the value i?c = /^^^^ ~ 
60 is not unreasonable. Clearly the calculation we have made is quite rough, but the 
exponent should not be very sensitive to details, in particular because k occurs only 

as its logarithm in (4). Notice also that the estimate k — K'J'^ is proposed instead of a 
fundamental calculation which is beyond current possibilities. Altogether, the agreement of 
(4) with experiment, with a plausible value of k, supports the physical picture of turbulence 
that we have presented. 
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